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ABSTRACT 

Microflares are small activities in solar low atmosphere, some are in the low corona, 
and others in the chromosphere. Observations show that some of the microflares are trig- 
gered by magnetic reconnection between emerging flux and a pre-existing background 
magnetic field. We perform 2.5D compressible resistive MHD simulations of magnetic 
reconnection with gravity considered. The background magnetic field is a canopy-type 
configuration which is rooted at the boundary of the solar supergranule. By changing 
the bottom boundary conditions in the simulation, new magnetic flux emerges up at the 
center of the supergranule and reconnects with the canopy-type magnetic field. We suc- 
cessfully simulate the coronal and chromospheric microflares, whose current sheets are 
located at the corona and the chromosphere, respectively. The microflare of coronal ori- 
gin has a bigger size and a higher temperature enhancement than that of chromospheric 
origin. In the microflares of coronal origin, we also found a hot jet (~1.8 x 10 6 K), which 
is probably related to the observational EUV/SXR jets, and a cold jet (~10 4 K), which 
is similar to the observational Ha/Ca surges, whereas there is only an Ha/Ca bright 
point in the microflares of chromospheric origin. The study of parameter dependence 
shows that the size and strength of the emerging magnetic flux are the key parame- 
ters which determine the height of the reconnection location, and further determine the 
different observational features of the microflares. 

Subject headings: Sun: microflare — Sun: magnetic reconnection — methods: numerical 



1. INTRODUCTION 



Microflares, or subflares, which are smal l-scale and short-lived solar activities, have already 

been studied for many years since last century ( Smith h Smithfl963 ; Svestkafl976 ; Tandberg-Hanssen h Emslie 
19881 ). The typical size, duration and total released energy are 5-20 arcsecs, 10-30 minutes and 
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2002; 
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2006, : 


2010), respectively. Microflares have been 


lany wavelengthes including Ha 


(Sc 


imieder et al. 


1997: 


Chae et al. 1999: Tane 


et al. 


2000). EUV ( 


iCmslie & Noves 


197J 




Porter et alj 


1984: 


Chae et al. 1999: 


Brosius & Holman 


2009; 


Chen & Dine 


2010), soft X-ray ( 


Golub et a . 


1974. 1977: iTane et al 


. 2000: Shimizu et al. 


2002: 


Kano et al. 201ol). hard X-rav ( 


Lin et al. 1984: 


Qiu et al. 


20041: lNinell2008: 


Brosius & Holman 


2009), 


and microwave (jGarv & Zirin 


1988: 


Gopalswamv et al. 


1994; ( 


3arv et al. 


1997). However, not all 



microflares have emissions at all wavelengths. Observations show that most of bright X-ray mi- 
croflares also appe ar at EUV and Ha bands, but only part of Ha microflares have their counterparts 
in X-ray emission (jZhang et al.ll2012l ). 



Observational characteristics of microflares, such as the heating, the relation with magnetic 
field, the duration and the coincidence between different wavelengths, etc, imply that mic roflares 
are produced by magnetic reconnection, similar to big flares. For example, iQiu et al.1 (J2004J) found 
that about 40% of microflares show hard X-ray emissi ons at over 1 keV and microwave emissions 
at about 10 GHz, typical features for flares. Recently, iNine) (|2008l ) found that roughly half of the 
microflares display the Neupert effect as previously revealed in f lares. On the othe r hand , some 
other microflares may have their origins in the lower atmosphere. iBrosius &; Holmanl (|2009i ) found 
that the microflares are bright in the chromospheric and transition region spe ctral lines , which 
is consistent with the chromospheric heating by nonthermal electron beams. iJess et al.l (|2010i ) 
concluded that the microflares in their study are due to magnetic reconnection at a height of 200 
km above the solar surface. 

Since a part of the microflares are located at emerging flux regions, they are probably due 
to magnetic rec onnection driven by th e new emerging magnetic flux (EMF) with a pre-existing 
magnetic field. ISchmieder et al.l (119971 ) found that the X-ray loops of microflares appear at the 
locations of EMF. 



Chae et al. 



(| 19991 ) found some repeatedly occurred EUV j ets where the pre - 
existing magnetic field was "canceled" by the new EMF with opposite polarity. iTang et al.l ([2000 ) 
indicated that the new EMF successfully emerged abo ut 20 minutes before the peaks of Ha and 
soft X-ray brightenings. Half of the events studied by IShimizu et al.1 (12002 ) s how the sma l l scale 
emergences of magnetic flux loops in the vicinity of the transient brightenings. iKano et al.1 (120101 ) 
found that EMFs and moving magnetic features (MMFs) are related to the energy release for 
at least half of the microflares around a well developed sunspot. Therefore, to understand the 
dynamics of microflares, it is crucial to simulate the magnetic reconnection associated with EMF. 



The numerical simulations of flux emerge nce scenario started from 1980s (jHoriuchi et al. 



Matsumoto et al 



expe r iments ( Maeara 2001 



1988: 



2006; Isobe et al 



Shibata et al 



200 



7: 



Fan 



1989a 



2001 



Hood et al 



1988 



bj). Later, more r e alistic and compl i cated numerical 



Manchester et al.l 12004 : llsobe et al.l 12003 : iLeake fc Arber 



2009) were presented to study the emergence of twisted flux 
from the convection region to solar corona. The simulations are mainly based on the Parker 



instability which is proposed by iParkerl (119661 ) and their results can be applied to account for the 
formation of filament or newly active region. More 2D and 3D numerical experiments are focused 
on the interaction between the emerging flux and the per-existing coronal magnetic configuration 
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to study the solar eruptions. With 2D simulations, it is found that the magnetic reconnection 
due to flux em e rgence can well explain solar jets or explosive events (jYokovama fc Shibatalll995l ; 
Jin et al.l Il996l ; iNishizuka et al. 2008 ; Ding et al.l I201U . l201ll ). or can even trigger the onset of 
coronal mass ejections (jChen &: Shibatal l200dl ) . Further 3D simulations gave the same but more 
detailed results, which revealed that the emerging flux plays a key role in the formation of many 
sol ar activities, such as type II spicules, erupting filaments or flux ropes, Ellerman bombs, and so 



on (IGa lsgaar d et al 



2005, 



Martmez-Svkora et al.ll201ll ). 



2007: Tdrdk et al 



2009 



Archontis Torok 



2008; 



Archontis fc Hood 



2009 



From the description of previous works, it is seen that EMF is responsible for many different 
phenomena, and one main cause of the divers ity of the dynamics is that the r econnection happens 



at diffferent heights in the solar atmosphere (jShibata 



1996 



Chen et al.l ll999). Our simulation in 



this paper is mainly focused on the microflares in a canopy-type magnetic configuration, where 



(Jiang et al. 


2010; 


Xu et al. 


2011 



we presented the 2.5 dimensional (2.5D) magnetohydrodynamic 
(MHD) simulation of magnetic reconnection in the chromosphere, which reproduces qualitatively 
the temperature enhancement observed in chromosphe r ic mi croflares. However, the magnetic con- 



figuration in that simulation, the same as in lChen et al.l (|200ll ). is too idealized and the reconnection 



site was specified at the computational center. Besides, no corona was included, so that the simu- 
lations can tell nothing about the EUV and SXR features. In this paper, we adopt a more realistic 
magnetic configuration and extend the solar atmosphere up to the corona in order to see what 
determines a microflare to be of coronal origin or of the chromospheric origin, and to compare the 
difference between the two types of microflares. The paper is organized as follows: the numerical 
method is described in Section [2j the numerical results are presented in Section [3l Discussion and 
summary are given in Section [H 



2. NUMBERICAL METHOD 



2.1. Basic Equations 



How to keep the magnetic field divergence free is one of the big problems in all MHD codes. 
Because of the discretization and nu merical errors, the performance of the MHD code can be 
unphysical (jBrackbill fc Barnesl I1980T). There are se veral ways to maintai n V • B = fo r MHD 



equations: fl) 8-wave fo rmulation (IPowell et al.lll999i), (2) the CT method (|Evans fc Hawlevlll988l : 
Stone Normanl Il992l ). (3) the projection scheme (jBrackbill &: Barnesl Il98dl ). The comparison 
betwe en these me thods showed that different methods have their own advantages and disadvan- 
tages (|T6thll2000i ). Besides, one can rewrite the original MHD equations by using vector potential 
A instead of the magnetic field B, or by using the vector magnetic potential or Euler potential. 
The advantage is that the divergence free condition is alw ays satisfied , howe ver, the MHD equa- 
tion should be rewritten. Our method, as introduced in Ijiang et al.l (|2012l ) adopts the 8-wave 
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Table 1: Normalization Units 



Variable 


Quantity- 


Unit 


Value 


T 


Temperature 


To 


10000 K 


P 


Density 


Po 


2.56 x 10~ 8 kg m~ 3 


x,y 


Length 


Lq = T K B /(mg s ) 


301.6 km 


P 


Pressure 


po = p ToK B /m 


2.12 N m" 2 


V 


Velocity 


vo = (po/po) 1/2 


9.09 km s _1 


B 


Magnetic field 


Bo = (poPo) 1/2 


16.3 G 


t 


Time 


to = Lq/vq 


33.1 s 



method. As suggested by iDedner et al.1 (120021 ) . we use the Extended Generalized Lagrange Mul- 
tiplier (EGLM)-MHD equations rather than pure MHD equations, which include two additional 
waves to transfer the numerical error of V • B. The local divergence error can be damped and 
passed out of the computational domain. The adopted dimensionless EGLM-MHD equations with 
resistivity and gravity included are given as follows: 



dp 

at 



+ V • (pv) = 



djpy) 
dt 



+ V • ( ( p + -B A ) I + /3vv - BB 



OB 



(V • B) B + pg 



— + V • (vB - Bv + VI) = - V x (r/V x B) , 
dt 



(1) 
(2) 
(3) 



de _ 



e + -B 2 +p) - B (B • v) 



-B-(V^)- V'((i)VxB) xB)+pg-v , (4) 



(5) 



where eight independent conserved variables are the density (p), momentum (pv x , pv y , pv z ), mag- 
netic field (B x , By, B z ), and total energy density (e). The expression of the total energy density is 
e = p/(7 — 1) + pv 2 /2 + B 2 /2 where 7 = 1.1 is taken in all our computations. The pressure p and 
temperature T are dependent on the eight conserved variables, g is the gravity vector, and rj the 
magnetic resistivity coefficient. Finally, a unity matrix I is the unity matrix. The main variables 
are normalized by the quantities given in Table HJ 



(Dedner et al. 


2002; 


Matsumoto 


2007) 



C h IS 



suggested by IDedner et al.l (|2002l ) , the expressions for Ch and c p are: 



- 5 - 



where At is the time step, Ax and Ay are the grid sizes, c c ji is a safety coefficient less than 1. 
c<j € (0, 1) is a problem dependent coefficient to decide the damping rate for the waves of divergence 
errors. We can see that Ch and c p are not independent of the grid resolution and the scheme used. 
Hence we have to adjust their values for different situations. 



2.2. Initial Condition 



The computational box is located in the x-y Cartesian plane. The x-axis is parallel to the 
solar surface while the y-axis is perpendicular to the photosphere. As shown in Figure [TJ the 
computational domain is —100 < x < 100 and < y < 200, where the length unit is Lq = 301.6 
km. In the simplified chromosphere and photosphere the temperature is set uniform with the value 
of 0.6, after being normalized by To = 10000 K. This layer extends from the bottom to y = 8.5. 
After a thin transition region with a thickness of 1, the plasma temperature rises rapidly to 100 in 
the corona, which corresponds to 1 MK. Given the temperature distribution we can get the density 
and pressure distributions according to the hydrostatic equilibrium. The initial distributions of 
density, gas pressure, plasma f3 (the ratio of gas to magnetic pressures) and temperature along the 
red line in Figure Q] are outlined by Figure [2j 

As shown in Figure [TJ two canopy-shaped magnetic structures are separated by a distance of 
100 (corresponding to 30 Mm) in order to mimic the t wo boundaries of a supergranu l ar cel l . The 
canopy magnetic field, which was originally proposed by iGabriell (|1976l ) and iGiovanellil (Il980l .ll982) 
for studying the chromospheric and photospheric magnetic fields, is rooted at the supergranular 
boundaries. In order to create a canopy magnetic field, we introduce several "magnetic charges" 
locating beneath the photosphere. Accroding to the Gauss's theorem in the magnetism, the mag- 
netic field (at the position r) generated by a "magnetic charges" (at the position r') can be written 
as B = ce(r — r')/|r 
of the field) and B = 



- r'| 2 in two dimensional geometry (c# is a constant related the strength 
cg(r — r')/|r — r'| 3 in three dimensional geometry, which is similar to the 



electric field. The locations of the "magnetic charges" are depicted in Figure [TJ The magnetic 
config urati on in our s i mulat i on is not exactly the same as the canopy model proposed by iGabriel 
(Il976h and iGiovanelhl (119801 . 119821 ). 
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rona 



Photosphere 

Fig. 1. — Initial condition for the magnetic configuration and temperature distribution. The ranges 
of the computational box is from —100 to 100 in the x-direction and from to 200 in the y-direction. 
The shadowy region means the damping zone. The "magnetic charges" are located beneath the solar 
surface. Note that the "magnetic charges" are samples to show the positions of these "charges" . 
The units of length, temperature, time and velocity are 301.6 km, 10000 K, 33.1 s and 9.09 km s _1 , 
respectively. 
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Fig. 2. — Initial distributions of the density (p), gas pressure (p), plasma beta (/3) and temperature 
(T). The density, gas pressure and plasma j3 distributions use the same left y-axis and the right 
blue y-axis is for the temperature. The units of density, gas pressure, length, temperature, time 
and velocity are 2.56 x 1(T 8 kg m" 3 , 2.12 N nr 2 , 301.6 km, 10000 K, 33.1 s and 9.09 km s _1 , 
respectively. 
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2.3. Resistivity Model 



In order to reproduce a fast reconnection (jPetschekl 1 1 964 ) . we adopt a non-uniform anomalous 
resistivity. This kind of resistivity may be due to some microscopic instabilities, although how 
these instabilities drive the macroscopic r econnection is still not clear. The mathematical form of 
the a ssumed anomalous resistivity model (jugail fl98~ 
200lh is: 



Chen fc Shibatall2n0f1 : lYokovama fc Shibata 



V 



mm 



IJI 

jc 



1 



, Vr, 



for \j\ > j c 



(8) 



where |j| is the total current density, j c the threshold, above which the anomalous resistivity is 
excited, and r\ max the maximum value of the resistivity. As shown by Equations ©, the resistivity 
(r/) is dependent on the current density. The anomalous resistivity (77) is switched off when |j| < j c . 



2.4. Boundary condition 



The magnetic flux emergence is realized by changing conditions at the lower boundary (jForbes &; Priest 



1984 



Chen &: Shibata! 12000 : iDing et al.ll2010l ). In our case, B x and B y , rather than the magnetic 



flux function in previous works, are independent variables, we change directly the magnetic field 
value at the boundary. The magnetic field (B x , B y , B z ) changes with time while the other variables, 
e.g., p, v, and p, are fixed to the initial values at this boundary. The mathematic form of the EMF 
is given by: 



B x ( X ,y) = -B e ±^2 , (9) 
B y (x,y) = B*^ZJ^ , (10) 

in the range x 2 + y 2 < r 2 , y < 0, where B e is the strength of the EMF, t e the end time of emergence, 
vq the half width of the EMF, x e and y e the center coordinates of the EMF in x- and y-direction. 
In our simulation esses, x c and y e are set to be and —3.6, respectively. Formulae (f9j)- (fT0|) are 
applied to the magnetic field B until t = t e . Later, all the physical values are fixed. 

Free boundaries (equivalent extrapolation) are adopted at the upper side y = 200. However, 
because of the low plasma beta (/3 = 2p/B 2 ) the upper boundary becomes unstable when the jets 
or waves propagate to this boundary. Thus we add a damping zone from y = 180 to y = 200 in 
y-direction as shown in Figure HJ The waves or jets will be damped to the initial value as the 
formulae given below: 



u(x, y) = u(x, y)D(y) + u (x, y) (1 - D(y)) 



(11) 
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where u represents varibles p, v, and p (not B), uq is the initial value, D the damping function, 
y s the start of position for the damping zone in y direction, y e the end position for damping zone. 
In all our cases, y s = 180 and y e = 200. Nonetheless, the damping zone may also lead to some 
reflections of waves or jets, but these reflections have little effect on the results since they are far 
away from the lower atmosphere which we are interested in. We use fixed boundary condition for 
the left and right boundaries at x = —100 and x = 100. 



2.5. Numerical Scheme and AMR Grid 



The code used in our simulation is MAP (j Jiang et al.ll2012l ). which is developed by the solar 
group of Nanjing University. The MAP code is a FORTRAN code for MHD c alculation with the 



adaptive mesh refinement (AMR) (jBerger fe Qligerill984j : iBerger &; Coleilalll989l ) and Message Pass 



ing Interface (MPI) parallelization. MAP have three optio nal numerical schemes f or the MHD part 



name ly, modified Mac Cormack Scheme (|Yu &: Liu 



2001 



Lax-Fridrichs scheme (Toth &; Odstrcil 



19961 ) and weighted essentially non-oscillatory (WENO) (IJiang fe; Wulll999l ) scheme. In this paper, 
we mainly use the WENO scheme, which has a higher accura cy than the other two. Moreover, it 
can keep the total variation diminishing (TVD) (|Hartenlll997l ) property for the MHD part without 
any additional artificial viscosity. The base resolution of our simulation is 256 x 256 and the mesh 
refinement level is 4 in all of our cases. Thus the effective resolution is 2048 x 2048 globally and 
the minimum grid size is 0.1 (about 30 km). 



3. RESULTS 

Our simulations show two typical cases, one with the reconnection occurring in the corona while 
the other mainly in the chromosphere and partly in transition region. The dynamics of magnetic 
reconnection processes is similar in the early stage in the two cases. Both show the following 
evolution: (1) an EMF emerges from the center of the supergranule; (2) the EMF reconnects with 
the pre-existed magnetic field; (3) temperature is enhanced near the X-point and the reconnection 
inflow and outflow appear. However, the final results are very different. The difference between 
these two cases are presented in 13.11 and 13.21 subsections. 



3.1. Coronal Microflare 



In this case, we set the parameters of the EMF to be B e = 32 (corresponding to 520 G), 
f]max =0.1 and ro = 8 (corresponding to ~ 2400 km). The reconnection occurs when the local 



- 10 - 



current density is larger than the threshold (J c = 10) in the formula ©. Figure [3] depicts the 
results of the simulation at different times. In this figure the color stands for the temperature, 
solid lines for magnetic field and arrows for velocity. The upper panel of Figure [3] shows that the 
EMF already rose into the photosphere and chromosphere. At the time 55, the transition region 
has been pushed up to the height y = 35 (i.e., about 10000 km), and fast reconnection has not 
happened yet. Later, when the magnetic reconnection rate reaches the maximum (time=70, as 
shown by the upper left panel of Figured]), a large amount of magnetic energy is released to the 
internal energy and kinetic energy. In this case, the reconnection occurs at the corona and the size 
of this microflare is about 50 from x = —40 to 10 (~ 15000 km, i.e. ~ 20 arcsec). The hot jet 
(~ 1.8 x 10 6 K) and cold jet (~ 10 4 K) ha ve for med around this time that are similar to the results 
simulated by lYokovama fc Shibatal (|1995l ) and iNishizuka et al.1 (|2008l ) . 



The reconnection rate and the one dimensional (ID) distributions of density, temperature 
and y-component of the velocity along the red line in Figure [3] are illustrated in Figure HI The 
reconnection rate is calculated by the ratio of the inflow speed (Vj n ) and the local Alfven speed 
(Va). Both two speeds are measured around the current sheet. The reconnection rate reaches the 
maximum around tim e=70 (38 minut es) and its value is about 0.09 that is within the value of 
0.01 - 0.1 predicted by iPetschek! (|1964l ). The high reconnection rate partly results fro m some small 



plasmoids ejections from the X-point, which may increase the reconnection rate (jShibata et al 



1995 



Shibatalll996l ). While after time=76, some numerical errors occur at the center of the EMF. 



Thus, the result after that time is unreliable although the reconnection rate increases again (see the 
upper left panel of Figured]). A roughly estimated life time of this coronal microflare is about 12 
minutes if we take the duration being from the beginning (time=55) to the ending (ti me=75) of the 



main reconnection. That is consist ant with the observations values (10 - 30 minutes: IShimizu et al 



2002: lFang et al.ll2006l . l2010l ). The other three panels show one dimensional distributions of density, 
temperature and y-component velocity at time=45, respectively. As we mentioned above, a hot jet 
and a cold surge are formed in this case. The hot jet (~ 1.8 x 10 6 K) which originates from the 
reconnection region, is ejected to the higher corona w ith the velocity of 16 (about 140 km s -1 ), 
which corresponds to the observational EUV/SXR jet ( Chae et al.1 19991 ; Brosius Holman 2009 ). 
The denser cold surge (~ 10 4 K and the density is two orders of magnitude larger than that of the 
hot jet), which is drawn to the corona by the hot jet. The cold surge falls to the lower atmosphere 
later with the spe ed of 3 (~ 30 km s^ 1 ). This process is similar to the observational Ha/Ca surge 
(jChae et al.lll999l ). The estimated size of the EUV/SXR or Ha bright point in our simulation is 
about 50 from x = —40 to 10 (~ 15000 km, i.e. ~ 20 arcsec). 



3.2. Chromospheric Microflare 

In this case, the strength and half width of the EMF are set to be B e = 12 (200 G), rj max = 0.1 
and ro = 6 (~ 1800 km). Compared the coronal case, the reconnection process is similar but 
the simulation result is totally differenct, since the coronal magnetic field prevents the weak EMF 
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Time=35.0 




-100 -50 50 100 

x 



Time=55.0 




-100 -50 50 100 

x 



T ime=70.0 




-100 -50 50 100 

x 



Fig. 3. — Temperature distributions (color scale), projected magnetic field (solid lines) and velocity 
field (vector arrows) at times 35 (upper panel), 55 (middle panel) and 70 (lower panel). The units 
of length, temperature, time and velocity are 301.6 km, 10000 K, 33.1 s, 9.09 km s -1 , respectively. 
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Reconnection Rate 



Time = 70.0 
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0.0100 
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0.0001 
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Time = 70.0 



Time = 70.0 
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-20 -10 10 20 30 40 



-20 -10 10 20 30 40 



Fig. 4. — The upper left panel shows the magnetic reconnection rate using the ratio between the 
inflow and Alfven speeds. Three vertical dashed lines in this panel outline the times used in Figure 
[3l The other three panels are one dimensional distributions of density (RO), temperature (TE) 
and y-component velocity (VY) along the solid red line shown in the lower panel of Figure El The 
horizontal dashed lines show the initial values while the vertical lines indicate the boundary between 
the hot and cold jets in these three panels. The density, length, temperature, time and velocity 
units are 2.56 x 1CT 8 kg m~ 3 , 301.6 km, 10000 K, 33.1 s and 9.09 km s" 1 , respectively. 
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emerging to the height as high as the coronal one. Therefore, the magnetic field accumulated 
in the chromosphere until the local current density exceeds the threshold (j c = 10, the same as 
the coronal case). Finally, we got the chromospheric microflare resulting from a fast reconnection 
mainly in the chromosphere and partly in the transition region. The result of this chromospheric 
microflare is given by three snapshots in Figure [5j Similar to Figure El the three panels correspond 
to the times before the reconnection, the time when reconnection rate is increasing, and the time 
when the reconnection rate reaches to the maximum (as shown by the upper left panel of Figure 
E]), respectively. In this case, the reconnection occurs at the chromosphere and the size of this 
microflare is about 15 from x = —10 to 5 (~ 4500 km, i.e. ~ 6 arcsec). There is no obvious 
reconnection jets ejected into the corona and we only obtain some temperature enhancements in 
the chromosphere and partly in the transition region. 

The reconnection rate and the one dimensional (ID) distribution of temperature along the red 
line in Figure [5] are illustrated in Figure [6j The physical process is much simpler than the corona 
case. As the reconnection occurs in the dense chromosphere, it is very hard to heat the plasma 
to a very high temperature. The right panel of Figure [6] is the temperature distribution at the 
the height of 1500 km. The inital value of the temperature is 0.64 (6400 K) and the temperature 
enhancement is about 800 K at this layer. The low temperature region is due to the expansion of 
the EMF. As we know, the plasma satisfies the frozen-in condition in the solar atmosphere without 
the resistivity. Therefore, the plasma of the lower atmosphere is dominanted by the strong magnetic 
field and finally such an adiabatic expansion will reduce the temperature. The low adiabatic index 
(7 = 1.1 in our simulations) can decrease such an effect, which also makes the simulated system 
more similar as an isothermal process. Although we can see two regions where the temperature is 
increased as seen in Figure (H they are not be observed as two Ha bright points since the hot region 
is a shell covering the EMF as seen from the lower panel of Figure [5j That is to say, we can only 
observe one Ha bright point if the line of sight is along the y-axis. The size of this bright point 
is about 15 from x = —10 to 5 (~ 4500 km, i.e. ~ 6 arcsec). There is no obvious temperature 
enhancement in the corona. Thus this case can explain the microflares with only Ha emssion. 
The inflow (downward) speed in the lower panel of Figure [5] is about 6 km s _1 which may be 
observed as a redshifted Doppler velocity in coronal lines. The outflow velocity of the reconnection 
flow in the chromosphere is about 20 km s _1 , but the direction of this outflow is redirected by 
background canopy-type field. Eventually, this coronal outflow (~ 4 km s _1 ) moves upward next 
to the downward inflow as shown by the lower panel of Figure [5j 



3.3. Parameter Dependence 

So far we have simulated one coronal case and one chromospheric case. An important question 
is which parameters determine whether reconnection happens in the corona or in the chromosphere 
when an EMF appears. In this subsection, we study three parameters, i.e., the strength of the 
EMF (B e ), the maximum value of resistivity (r] max ), and the half width of EMF (ro). In order to 
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Time=40.0 
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Time=72.0 
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Fig. 5. — Temperature distributions (color scale), projected magnetic field (solid lines) and velocity 
field (vector arrows) at time 20 (upper panel) , 40 (middle panel) and 72 (lower panel) . The units of 
length, temperature, time and velocity are 301.6 km, 10000 K, 33.1 s and 9.09 km s _1 , respectively. 
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Reconnection Rate Time = 72 




Fig. 6. — The left panel shows the magnetic reconnection rate using the ratio between the inflow 
and Alfven speed. Three vertical dashed lines in this panel outline the times used in Figure [5j The 
other panel are one dimensional distributions of temperature (TE) along the solid red line shown 
in the lower panel of Figure [5j The horizontal dashed lines show the initial value. The length, 
temperature and time units are 301.6 km, 10000 K and 33.1 s, respectively. 



understand the effects of different parameters, we perform extensive simulations by changing one 
parameter with others being fixed. 

Figure [7] shows the dependence of the reconnection process on different strengths of the EMF 
magnetic field. Four values of B e are studied, i.e., B e = 12,22,32 and 42, whereas the other two 
parameters are fixed, i.e., rj max = 0.1 and tq = 8. With the increasing of the strength of the EMF, 
the magnetic reconnection rate can reach the maximum much faster and the maximum reconnection 
rate becomes larger and larger when we set a stronger initial magnetic field. The right panel shows 
how the maximum temperature and the height of the current sheet center change with the strength 
of EMF magnetic field. The Max Temperature and Height is very sensitive to the parameter B e . A 
stronger EMF can result to a higher altitude which can generate a stronger current density. Thus, 
it is easy to trigger the reconnection at an earlier time. The reconnection is faster and the outflow 
is hotter than a weak EMF. When we set B e = 12, we find that the height is 12.5 (3750 km), which 
is located nearly above the transition region. If we keep reducing the value B e , the reconnection 
will occur in the chromosphere. 

Figure [8] shows the dependence of various quantities on different values of maximum resistivity 
(rjmax)- Four different values are selected, i.e., rj max = 0.01,0.05,0.1 and 1.0. The other two 
parameters are fixed, i.e. B e = 32 and ro = 8. We got the similar results as described in our 
previous paper (| Jiang et al.ll2010i ). The variation of the resistivity value does not change the results 
too much. The reconnections almost occur around the heigh of 28 (8500 km) and the temperature 
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Fig. 7. — The left panel shows the dependence of Max Vi n /VA (maximum magnetic reconnection 
rate) and Max Rate Time (the time when the reconnection rate reaches the maximum) on different 
strength of the EMF magnetic field. The right panel depicts the dependence of Max Temperature 
(the maximum temperature at the Max Rate Time ) and Height (the height of current sheet center 
at the Max Rate Time ) on different strengths of EMF magnetic field. The length, temperature, 
time and magnetic field units are 301.6 km, 10000 K, 33.1 s and 16.3 G, respectively. 




Fig. 8. — The left panel shows the dependence of Max Vi n /VA and Max Rate Time on different 
maximum values of resistivity (rj max ). The right panel depicts the dependence of Max Temperature 
and Height on different maximum values of resistivity (rj max ). The length, temperature, time and 
resistivity units are 301.6 km, 10000 K, 33.1 s and 3445.1 SI m, respectively. 
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enhancement at such a height is about 350 (3.5 x 10 6 ). 
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Fig. 9. — The left panel shows the dependence of Max Vi n /VA and Max Rate Time on different half 
width (ro) of the EMF. The right panel depicts the dependence of Max Temperature and Height 
on different half width (ro) of the EMF. The length, temperature, time and length units are 301.6 
km, 10000 K, 33.1 s and 301.6 km, respectively. 

The dependence on different half widths (ro = 4, 6, 8 and 12) are depicted in Figure [H The 
other two fixed parameters are B e = 32 and r\ max = 0.1. From the reconnection start time to the 
maximum time, a bigger EMF can rise to a higher altitude. We know that in the higher altitude 
the plasma f3 is lower, thus the magnetic reconnection can be very fast with the same resistivity. 
At the same time, the strength of the magnetic field of EMF becomes weaker and weaker after 
expanding to such a long distance, so the temperature enhancement becomes smaller and smaller 
as indicated by Figure Therefore, the left panel of Figure [9] shows the reconnection rate and 
the reconnection duration increase with the increasing half width of the EMF. In the right panel, 
we find that the height of X-point is also sensitive to the size of the EMF. Combining the small 
magnetic strength of the EMF, we can get the result of chromospheric microflare as described in 
Section I 



4. DISCUSSION AND SUMMARY 

In this paper, we simulate the microflares produced by EMFs. The EMFs with different 
strength and size lead to different results. Two types of the microflares in our simulations, i.e., of 
coronal origin and of chromospheric origin, can be illustrated by the cartoon in Figure [TUl According 
to our simulations, small microflares with weak Doppler velocities and only Ha emissions are most 
likely due to the reconnection in the lower solar atmosphere (for instance, the chromosphere or 
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photosphere). However, big microflares with obvious Doppler velocities and emissions in both low 
and high temperature wavelengthes (Ha and EUV/SX R) originate fr o m the reconne c tion i n the 
corona. Th i s mod el is similar to the unified model of IShibatal (| 19961 ) ; IShibata et al.l (|2007l ) and 
Chen et al.l (|1999l ) in the sense that the different behaviors of eruptions are determined by the 
different height of the magnetic reconnection. According to the result in this paper, the strength 
and the size of the EMF determines whether the microflare is of coronal origin or of chromospheric 
origin when the EMF emerges from the center of a supergranule, probably dragged by the convective 
upflow. Of course, EMFs may also appear at other sites of a supergranule, which would also affect 
the height of the reconnection site. 
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Fig. 10. — The cartoon shows the basic physical process of our simulated two types of microflares, 
i.e. of chromospheric and coronal origins. 

As mentioned in Section 13.11 we find some plasmoids in our simulations. Figure [11] shows 
the current density distribution of a small region [-35, 0] x [10, 45] in Figure [3j From both the 
temperature and current density distributions (Figure 131 and [TTj) , we can see a plasmoid is formed at 
the center of the c urrent sheet. The ejection of the plasmoid can incre ase the magnetic reconnection 
rate temporarily ( Shibata et al. 1995 ; Shibata 1996 : Shen et allboill ). In the reconnection process, 
the plasmoids are generated one by one, therefore we can see some oscillations of the reconnection 
rate in Figure HI However, for the chromospheric case, the Alfven speed is much smaller than 
the coronal o ne, which leads to a smaller magnetic Reynolds number. Thus, similar to the results 
presented by IShen et al.l (|201ll ). there is no plasmoid produced in the chromospheric case. The 
corresponding reconnection rate in Figure [6] also shows less oscillations. In the coronal case, as 
the reconnection occurs in the corona, both hot jets and cold surges are ejected. No doubt that 
the EUV/SXR emissions from the hot jet will be observed first, whereas Ha/Ca bright points and 
surges appear late r. If we take the tem perature response at the height of 1500 km, which is the Ha 
formation height (jVernazza et al.lll98ll ). the roughly calculated time delay between th e EUV/SXR 



and Ha emissions is about 3-5 minutes, which is comparable to the observations (jZhang et al 



2012J). We have to mention that the temperature enhancement and the time delay can be more 
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realistic if we include the thermal conduction, radiation, the high energe particles (non-thermal 
electrons and ions), and the effect of partially ionised plasma. For simplicity and for the sake to 
see how emerging flux with different sizes and strengths influencesthe reconnection, we kept the 
threshold of the current density j c in the anomalous resisitivity in Eq. (J8|) to be constant. In reality, 
the threshold of j c in the chromosphere might be higher because of its high plasma density. In order 
to see how j c would change the reconnection process, we performed other numerical experiments 
with higher j c in the chromosphere, and it is found that the basic results are the same, the main 
difference is that the commencement of reconnection is postponed only for the case th at the current 
sheet located at the chromosphere, as demonstrated bv lYokovama fc Shibatal (| 19941 ). 
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Fig. 11. — A amplified view of z-component current density (JZ) distribution of the conoral case at 
the time 70. The overall view of temperature distribution is shown by the bottom panel of Figure 
[3] at the same time. The range of this box is [-35, 0] x [10, 45], in which a plasmoid in the center 
of the current sheet is clearly seen. 



In summary, we successfully simulated the coronal and chromospheric microflares by using 
MHD simulations with a canopy-type magnteic configuration. The different observational behav- 
iors between coronal microflares and chromospheric ones are due to the height of magnetic recon- 
nection, which is determined by the size of the emerging flux, i.e., smaller emerging flux produces 
chromospheric microflares and larger emerging flux produces coronal microflares. Correspondingly, 
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the resulting microflares have different sizes. In the case of coronal origin, the sizes of the simulated 
microflares are 11000 — 16000 km, i.e., 15 — 22", which correspond to big microflares. We find a hot 
jet (~1.8 x 10 6 K for the typical case), which should be relates to the observational EUV/SXR jet, 
and a cold jet (~10 4 K for the typical case), which corresponds to the observational Ha/Ca surge 
or brightening. Some plasmoids generated at the center of the current sheet are ejected, which 
can increase the magnetic reconnection rate. In the case of chromospheric origin, the sizes of the 
microflares are 4200 — 4500 km, i.e., ~6", which correspond to relatively small microflares. As the 
reconnection occurs in the chromosphere, only the Ha/Ca brightenings show up in this case, where 
no significant SXR brightening and no plasmoids are produced in this case. The parameter survey 
qualitatively shows that the size and strength of the EMF are the key parameters which determine 
the height of the reconnection X-point. For some typical values, we can get the reconnection either 
in the corona or in the chromosphere. The parameter r} max has little effect on the final results. 

This paper deals with the microflares related to EMFs. It should be mentio ned that some mi- 
croflares are due to other mechanisms, e.g., the reconnec t ion b etween the MMFs (jHarvey Sz Harvey 



19731 ) and the pre-existing ma g netic field. iPriest et al.l (119941 ) proposed a canceling magnetic fea- 
tures model. iBrowning et al.l ((2003) performed 3D MHD simulations to study the role of kink 
instability in the magnetic energy releas process, which may also lead to microflares. It is also 
noticed that the simulated emerging flux in this paper does not includ e the twisted fi e ld, which is 



required for the emerging flux to survive the subsurface convection (s ee iBabcock 



1961 



1975 



FanJ|2009j, and references therein), e.g., iMurrav &: Hoodl (J2008) found that the tension force 



Piddington 



of twisted tube plays a key role in determining whether the flux can successfully emerge to the 
upper solar atmosphere in 3D numerical experiments. In the 2D coronal simulations, the twisted 
field can be easily considered by including B z , which would not change the results greatly. More- 
over, a single-fluid model for the solar atmosphere was adopted in this paper for simplicity. In 



reality, the chromosphere, especially the lower part, is weakly-ionzed, and pa rtial ionization s 



lould 



be considered, which wo u ld lead to strong ani sotropic Cowling resist ivity (ICowling. T. G 



Khodachenko et al 



Soler et al.ll2009l : ISingh et al 



2004 



Leake &: Arber 



2011 



2006), ambipolar diffusion (jVishniac &: Lazarian 



1957; 



1999; 



Moreover, the ionization of neutral atoms would consume a 



signi ficant part of the released energy in the magnetic reconnection (jChen et al.ll200ll ; iJiang et al 
2O10h . 
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